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Введение. В работе описаны уточненные расчетные модели 
трибосистем с учетом зависимости вязкости и проницаемости 
пористого покрытия от давления при частичном заполнении 
рабочего зазора 
расчетные модели позволяют выполнить сравнительный анализ 


смазочным материалом. Полученные 
вновь полученных и уже имеющихся результатов, что 
подтверждает большую приближенность новых моделей к 
реальной. 

Материалы и методы. На основе уравнений движения жидкого 
смазочного материала для случая «тонкого слоя» и уравнения 
Дарси определены значения поля скоростей и давления в 
смазочном слое. В пористом покрытии на поверхности шейки 
вала определены основные рабочие — характеристики 
подшипников скольжения. 

Результаты исследования. Разработана математическая модель 
радиального подшипника, позволяющая на основе численного 
анализа определять поля скоростей, несущую 


способность и силу трения с учетом зависимости от ряда 


давление, 


дополнительных факторов. 

Обсуждение и заключения. Полученные результаты могут 
быть использованы в учебном процессе и в инженерной 
практике трибосистем машиностроения. 


Ключевые слова: радиальный подшипник, гидродинамика, 
вязкость, проницаемость, неполное заполнение. 
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Введение. В настоящее время развитие машиностроения характеризуется увеличением в современных 


двигателях мощности поршневых и роторных машин. Повышение требований к их надежности и долговечности 


приводит к неизбежному росту нагруженности и сопряжений трибосистем. Ключ к решению этой проблемы лежит в 


совершенствовании конструкций узлов трения трибосистем. Наличие смазочного материала в трибосистемах 


обеспечивает режим жидкостного трения, при котором потери достаточно малы, а износ трибоузлов минимален. При 


разработке математических расчетных моделей и конструировании новых, а также при модернизации разработанных 


трибосистем необходимо уделять повышенное внимание происходящим процессам и явлениям, к которым относятся 
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механические, тепловые и химические процессы. Возникает необходимость в разработке новых расчетных моделей 
трибоузлов, обладающих большой маслоемкостью и демпфирующими свойствами [1-13]. 


Цель работы. Целью настоящей статьи является разработка расчетной модели трибоузлов с пористым 
покрытием при частичном заполнении рабочего зазора с учетом зависимости проницаемости и вязкости от давления. 


Постановка задачи. Рассматривается установившееся движение трибосистемы с пористым покрытием. Вал 
вращается с угловой скоростью (©, а подшипник неподвижен. Пространство между валом и подшипником не 
полностью заполнено смазочным материалом (рис. 1). 








Рис. 1. Рабочая схема 


Ето. 1. Кипспопа[ Фаэтат 


В полярной системе координат с полюсом в центре вала уравнение контуров вала с пористым покрытием С, 
и С иподшипниковой втулки С, задается в виде: 
® приий тт. ® ! ий ® 
Су: г=п-Н; С: Г=К; 


С: "=н.(1+Н)=к +есоз6; 


Неон 26 а (1) 
2 7 
где и — радиус подшипника скольжения; 7, — радиус вала с пористым покрытием; Н — толщина пористого слоя; 
Н — толщина смазочного материала; е — эксцентриситет. 
Зависимость вязкости и проницаемости задается в виде: 

и’ = ше”,  К=Ке®, (2) 
где и,— характерная вязкость, А, — характерная проницаемость пористого слоя; а — экспериментальная 
постоянная величина, р’ — гидродинамическое давление в смазочном слое, и’ — коэффициент динамической 
вязкости смазочного материала, К’ — проницаемость пористого покрытия. 


Исходные уравнения и граничные условия. 
Исходными уравнениями являются уравнения движения вязкой несжимаемой жидкости для случая «тонкого 
слоя», неразрывности и Дарси с учетом (2). 








2 ! ! ! 
а 
ду и 00 ог 00 
д ! | д / | д , 
5 т Ш 
ОГ и ФУ Г’ 00 


Связь между размерными и безразмерными переменными задается в следующем виде: 


(3) 








— В смазочном слое: 
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0 ! ! ' * ' 
и=—, и=Оци, у=Обу, р=рр, и=ши, 
Р 
(4) 
Ц Ох? 
, * 0-0 
РАНО, АА, РВ 
о 
— в пористом слое: 
ДН па Ро Р, (5) 
Здесь и’,у’ — компоненты вектора скорости, р’ — гидродинамическое давление в смазочном слое, Р’— 


гидродинамическое давление в пористом слое. 
С учетом (4) и (5) уравнение (3) примет следующий вид: 
би _ 1 фр 9 ди, ОР 1Р 1ОР_ 





















































Е —+—=0, — ——= 6 
д" е? 00° г 00 д" м д м” 00° (о) 
И соответсвующими граничными условиями: 
Ор 
и А = =— ‚ И=О0, у=0 при г=1[+1с050, 
| г=0 ди" № 
БО- (7) 
(6, )= р(0.)=0, р=Р *_70 > я = 0, 
Н ОР [И 
Н 
где М = ШИ = 
НИ 
Решение функции тока (6), с учетом граничных условий (7), будем искать в виде: 
д д о = 
и= а +0(,-,0), ›= — +7(+-, 0), 0(-.0)=#(), (+, 0) =-$ (©) 
р 
Г а ар {0 С Со 
Е, и(г, 0) = $(8), —--=е” т (8) 
1(0) (,6)- 9) 40 и? (0) #“(6) 
Подставляя (8) в (6), с учетом (7), получим: 
4 __ а41_. % са Г. - 
“=&, 2 — с, —46—— =0, (ЕЕ = $(0) (9) 
4 4 4 4 0 
гы а д ы = 
1(0)=1, 70)=-м- ‚ р(0,)=р(,)=0, #(@)=0, 7()=0. (10) 
Р о 
Решение задачи (9)—(10) находится интегрированием. В результате: 
ду _ 5. (2 а 
м Е ие) =а —-|— +1 Е +1. п 
ие) ЧЕ (1) 
Определение гидродинамического давления. 
Для определения безразмерного гидродинамического давления в смазочном слое имеем: 
о 2 С. 
ФЕ, 12 
4 1?(0) 13(0) > 
где #(0) =1+1с0$6. 
Введем обозначение: 
2" (13) 
и, продеффиринцировав обе части равенства по 0, получим: 
а —ме “Р ре —ме “Р ы Е (14) 
49 49 40 49 
С учетом (14) уравнение (12) запишется в виде: 
ЧЕ С с) 
— = -—а + (15) 
4 [55 ны 
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Интегрируя уравнение (15) с точностью до членов второго порядка малости О(т?), с учетом 


2(0,)=2(60,) =1, получим: 


2 =-о( (90-2 зш 0)-+ 5, (0 — Зизт 0))-+1 (16) 
Используя граничные условия 2 (0,) = 2(60, ) =1 будем иметь: 
с» =—С1| 1+ ИРИ (17) 
9—9, 


С учетом (17) уравнение (16) примет следующий вид: 


2=—@С; РИН +1 (18) 
9, —0, 
ИЛИ 
пост 
е “Р =-аст зт0- ны 19 
п 0,6, (19) 


Выполняя аналитическое разложение функции ев ряд Тейлора с точностью до членов 0(о—) 


включительно, получим: 
20со5 ый - ы ут 9. - ы 
=сп| 90 -—“ 20 
Р=ст 0,6, (20) 





С учетом (20) давление фильтрующегося смазочного материала в пористом слое будем искать в виде: 


20соз р 9 р 


р(",0)= ("Ем ре (21) 
- 1 





Подставляя (21) в уравнение Дарси (6) для определения в(” ) ‚ приходим к следующему дифференциальному 








уравнению: 
д) В В 
В"(к ) + : = =0 (22) 
Г 
С граничными условиями: 

0 А сз) 

4’ В Н 
Н 


Непосредственное интегрирование (22) с учетом (23) для функции К (") позволяет получить выражение: 


; Ни” о-2Ни +Н?* 
В(, ВАА ЕКА 2 (24) 
2% -2Н-+Н“ Н|\?2к-2Нь+Н7у 


Интегрируя уравнение неразрывности по & от 0 до 1, приходим к следующему уравнению: 





1 
=. (ЕЕ (25) 


С учетом (11), (21) и (24) для с получим следующее выражение: 


с = вл-амн(1- — ] (26) 





чи 0 
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р а "9 


2 
где А 
9, — 0, 


Тогда, с учетом с, окончательно для гидродинамического давления получим: 








202 
р= ы| —12 ВМ {шт 0- 40) + тв _ = | (27) 


Определение несущей способности и силы трения. 
С учетом (11) и (28) для составляющей вектора поддерживающей силы и силы трения получим выражения: 


К =Р%\ Рут 940 = 6' вов. А, —0. +2з1(0, +60, )зв1(6, — 0, 
0, 
-_ 2.4(0, + 0 т аи м ь. 2А И (28) 
й 2 2 _ 
2 
гов 240. +0 1) 9 5 ый мо о (9, _ 6, у 


К =Р%\ `1Рсоз04 = бт ее (в, + 0.) (©, —60,)- 
0, 


—А р о о + 
2 2 2 2 


Говми (© _0 р с0$ е —— ый УИ 9, р —2 до 96] 














"5 6, 0, —09, +4 т ——— 


(8) #6} 


0> [2 —у 2 -& (1+ Аи) 
Ей и | 1% (6) а А и еореб 0, +6, ыы и * 
2 й, 





И’ 1 
0,23 
0,22 
0,21 
0,20 
0,19 
а 
0 0,2 0,4 0,6 0,8 1 
Рис. 2. Зависимость несущей способности от параметра а, Рис. 3. Зависимость несущей способности от параметра а, 
характеризующего зависимость вязкости от давления и характеризующего зависимость вязкости от давления: 
параметра Л, характеризующего проницаемость [—М= 0,25;2 — М=0,5 


пористого слоя " | 
Ето. 3. Веатиз геязяапсе — рагатеег о геаНоп 


Ето. 2. Веатие геязяапсе — рагатеег о геаНоп спагачетлия Ше у15со5Йу — ргеззиге гёйаНоп: 
спагачетилия йе у15со5Йу — ргебзиге апа рагатаег М [— М = 0,25; 2 — М = 0,5 
геайоп, спатаметлив рогоиб [ауег регтеаБИйу 
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Выводы: 

1. С увеличением параметра а, характеризующего зависимость вязкости от давления, несущая способность 
подшипника имеет экстремум, что приводит в этой области к снижению силы трения. 

2. С увеличением параметра №, характеризующего проницаемость пористого слоя, несущая способность и 
сила трения резко уменьшаются. 
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